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Abstract 

Density  functional  calculations  have  been  earned  out  for  twelve  molecules  by  several 
different  Kohn-Sham  procedures  in  order  to  test  the  effectiveness  of  two  approximate  molecular 
energy  formulas  that  have  been  applied  earlier  in  the  context  of  Hartree-Fock  theory.  The  equation 

Emoiecular  =  |  (Vne  +  2Vnn) ,  in  which  Vne  and  Vnn  are  the  nuclear-electronic  attraction  and 

nuclear  repulsion,  reproduces  Kohn-Sham  energies  with  an  average  deviation  of  less  than  0.70%. 
Emoiecular  =  n j  ,  where  n;  and  are  the  molecular  orbital  occupancies  and  energies, 

gives  an  average  deviation  of  1.6%,  with  a  determined  empirically  for  each  functional 
combination.  The  relationship  3Emolecule  =  1.14]>>i  ^  +  2Vm  +  Vne,  an  analogue  of  an  exact 

Hartree-Fock  expression,  predicts  Becke-3  Kohn-Sham  molecular  energies  to  within  an  average 
0.23%  of  the  actual  values. 


1.  Background 

The  Hellmann-Feynman  theorem  can  be  written  in  the  form  [1-3], 


3E 


8Z 


A  7n 


=  V0, A 


(1) 


in  which  E  is  the  total  energy  of  an  N-electron  atom  or  molecule  and  Vo, a  is  the  electrostatic 
potential  at  nucleus  A  due  to  the  electrons  and  other  nuclei: 


V 


0,A 


■2 

B*A 


Zb 

Rb-Ra 


r  P(r)dr 

|r-RA 


(2) 


Zb  is  the  charge  on  nucleus  B,  located  at  Rb,  and  p(r)  is  the  electronic  density. 


Eq.  (1)  suggests  that  it  should  be  possible  to  express  energies  rigorously  in  terms  of 
electrostatic  potentials  at  nuclei,  and  such  formulas  have  indeed  been  derived  [4]. 
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The  molecular  formula,  eq.  (4),  has  the  notable  feature  that  it  is  simply  a  summation  over  atomic 
contributions,  with  no  explicitly  interaction  or  mixing  terms. 

The  quantity  ZVo  in  eq.  (3)  is  the  nuclear-electronic  attraction  energy,  Vne,  while 
XZAv0,A  in  eq.  (4)  equals  Vne  +  2Vnn,  Vnn  being  the  nuclear  repulsion.  It  follows  from  the 
A 

virial  theorem  that  the  electronic  repulsion  energies  are  included  in  the  integral  terms  in  eqs.  (3)  and 
(4).  Thus,  these  two-electron  quantities  can  be  evaluated  rigorously  from  the  one-electron 
properties  Vo  and  Vo,  a- 

Within  the  framework  of  the  Thomas-Fermi  model  [5-7],  which  is  an  early  and 
approximate  version  of  density  functional  theory, 

E„tom  =  |  ZV0  =  |  Vne  <5) 


Eq.  (5)  can  be  regarded  as  a  simplified  form  of  eq.  (3).  Using  Hartree-Fock  values  for  Vo,  eq.  (5) 
reproduces  the  Hartree-Fock  atomic  energies  to  within  an  average  of  about  2%  [4,8].  In  the  spirit 
of  eq.  (4),  the  extension  of  eq.  (5)  to  molecules  by  just  summing  over  the  atomic  contributions  was 

tested  [9]: 


Emolecule  =  rIZAVOtA  “  “(Vne  +  2Vnn) 

'  A 

The  average  error  obtained  with  eq.  (6),  relative  to  the  Hartree-Fock  energy,  is  in  the 
neighborhood  of  1%  [9]. 

In  Hartree-Fock  theory,  the  energy  of  a  molecule  can  be  expressed  exactly  by  [10], 


Emolecule  —  "**  ^nn  ^ee 


(7) 


where  the  nj  and  <E[  are  the  occupancies  and  energies  of  the  molecular  orbitals.  Eq.  (7)  can  be 
combined  with  the  virial  theorem  to  give, 


3Emolecule  “  Xniei  +  Vne  +  2Vnn  ^ 

i 

If  eq.  (6)  is  inserted  into  eq.  (8),  the  result,  as  shown  initially  by  Ruedenberg  [11],  is, 

E molecule  =  n  iei  ^ 

i 

Eq.  (9)  would  be  exact  within  the  framework  of  Hartree-Fock  theory  if  eq.  (6)  were  rigorously 
correct.  In  practice,  Ruedenberg  found  eq.  (9)  to  be  satisfied  to  within  2  -  4%  [11]. 

While  the  accuracies  obtained  with  eqs.  (6)  and  (9)  are  clearly  not  sufficient  to  permit  the 
direct  calculation  of  reaction  energetics,  these  relationships  have  aroused  considerable  interest,  on 
both  conceptual  and  applied  levels  [12-46].  Accordingly,  we  have  now  examined  the  validity  of 
eq.  (6)  in  the  context  of  Kohn-Sham  density  functional  methodology.  Since  eq.  (7)  is  no  longer 
applicable,  neither  is  the  derivation  of  eq.  (9);  however,  we  have  investigated  empirically  the 
possibility  of  replacing  the  factor  ^  *n  ecf  (9)  by  some  other  constant. 


2.  Procedure 

We  have  carried  out  Kohn-Sham  calculations  for  12  molecules  of  various  types,  using  five 
different  combinations  of  functionals  and  basis  sets.  The  former  include  the  Becke  (B)  [47],  the 
Becke-3  (B3)  [48],  the  Perdew-86  (P86)  [49],  the  Lee/Yang/Parr  (LYP)  [50]  and  the 
Perdew/Wang-9 1  (PW91)  [51].  The  basis  sets  were  the  6-3 1+G**  and  6-311++G(3df,3pd).  For 
purposes  of  comparison,  we  also  include  Hartree-Fock  treatments  at  two  levels,  HF/3-21G  and 
HF/6-3 1+G**.  All  geometries  were  optimized  by  the  respective  procedures,  using  the  Gaussian 
94  code  [52].  Finally,  in  order  to  verify  that  our  findings  extend  to  larger  molecules  as  well, 
HF/6-3 1+G**  and  B3P86/6-31+G**  calculations  were  carried  out  for  benzene,  aniline  and 
nitrobenzene. 


3.  Results 

Table  1  lists  the  key  computed  energy  quantities  obtained  for  twelve  molecules  by  seven 
different  procedures.  The  extent  to  which  eq.  (6)  is  obeyed  in  each  instance  is  shown  in  Table  2, 
while  the  ratios  Emolecule/2ni  ei  * which  are  relevant  to  eq.  (9),  are  in  Table  3.  The  results  for 

benzene,  aniline  and  nitrobenzene  are  given  in  Table  4. 


It  is  seen  in  Table  1  that  eq.  (6)  generally  overestimates  the  magnitudes  of  the  Hartree-Fock 
total  energies,  whereas  the  Kohn-Sham  are  overestimated  and  underestimated  to  roughly  the  same 
extent.  However  Table  2  shows  that  the  average  deviations  are  less  than  1%  of  the  correct  values. 

It  is  notable  that  they  are  little  affected  by  increasing  the  basis  set  or  by  changing  the  functional 
combination. 

If  eq.  (9)  were  satisfied  exactly  by  the  Hartree-Fock  data,  the  ratio  Emoiecuje/^nj 
would  be  1.500.  Table  3  shows  that  it  actually  ranges  between  1.48  and  1.64;  the  average  is  1.55, 
as  was  pointed  out  earlier  by  Ruedenberg  [1 1].  For  the  Kohn-Sham  results,  it  can  be  anticipated 
that  this  ratio  will  be  greater  than  for  the  Hartree-Fock,  since  the  Kohn-Sham  orbital  energies  are 
known  to  be  considerably  smaller  in  magnitude  [53].  Indeed  the  BLYP  values  of 
EmolecuIe/Xni  ei  3X0  between  1.75  and  1.85;  for  the  B3  functional  combinations  they  are 

somewhat  lower,  between  1.69  and  1.80,  perhaps  reflecting  the  partial  inclusion  of  the  Hartree- 
Fock  exchange  energy  [48,52]. 

Table  4  confirms  that  benzene,  aniline  and  nitrobenzene  fit  in  well  with  the  smaller 
molecules,  in  terms  of  both  eq.  (6)  and  the  Emolecule/^nj  ratio.  The  degree  to  which  the 

Hartree-Fock  data  satisfy  eqs.  (6)  and  (9)  is  particularly  striking. 

4.  Discussion 

As  was  pointed  out  earlier,  if  eq.  (6)  were  satisfied  exactly  in  Hartree-Fock  theory,  then  eq. 
(9)  would  be  as  well.  For  example,  the  HF/3-21G  and  HF/6-31+G**  results  for  diazirine  deviate 
from  eq.  (6)  by  only  0.13%  and  0.19%  (Table  2),  and  the  ratios  Emoiecuie/^nj  ej  are 
correspondingly  very  close  to  the  theoretical  1.500,  at  1.512  and  1.514  (Table  3).  In  fact  the 
Hartree-Fock  data  show  excellent  linear  correlations  between  the  percent  deviation  from  eq.  (6)  and 
the  ratio  Emolecule/Xni  ei  (Figure  1);  as  the  deviation  approaches  zero,  the  ratio  approaches 

1.500.  The  least  squares  plots  for  the  3-21G  and  6-3 1+G**  basis  sets  can  almost  be  superposed, 
and  they  predict  1.505  and  1.504,  respectively,  as  the  values  for  the  ratio  at  zero  deviation. 

We  have  shown  empirically  that  an  equation  of  the  form, 

Emolecule  =<*^ni  ei 

i 

is  a  reasonable  approximation  in  Kohn-Sham  theory,  accurate  to  within  about  1.6%.  At  the  BLYP 
level,  a  is  about  1.79,  and  for  the  B3  functional  combinations  it  is  about  1.73  (Table  3).  (The  four 
B3  values  for  each  molecule  are  nearly  the  same.)  However  in  the  absence  of  an  equivalent  of  eq. 
(7)  or  eq.  (8),  there  appears  to  be  no  reason  to  anticipate  a  correlation  between  the  deviation  from 
eq.  (6)  and  the  ratio  Emolecule  /£ni  ej .  It  is  therefore  particularly  interesting  that  such 


correlations  do  exist  (Figures  2  and  3).  The  five  plots  are  quite  similar  in  appearance,  although  the 
BLYP  is  somewhat  displaced  from  the  others  in  the  vertical  direction. 

The  existence  of  these  relationships  suggests  that  attention  be  focused  upon  the  value  of 
Emolecule  /Xni  €i  ^at  corresPonc*s  t0  zer0  deviation  from  eq.  (6).  For  the  B3  results,  this  is 

about  1.717;  for  the  BLYP,  it  is  1.772.  Since  eq.  (6)  is  satisfied  at  these  points,  the  derivation  of 
eq.  (9)  can  be  reversed,  using  the  appropriate  factor  instead  of  %,  to  arrive  at  analogues  of  eqs. 

(7)  and  (8).  For  example,  for  the  B3  functional  combinations, 

Emolecule  =  s,  =  f(1.145)5>i  *i 

i  Z  i 

■jEmolecule  =  l-145^nj  ej 


“Emolecule  +  TEmolecule  “  1.145^ nj  Gj  +  Emolecule  L145]Lni  6i  +  2Vnn  +  Vne 


3Emoiecule  “  l-145^ni  Gj  +  2Vnn  +  Vne 


(ID 


Emolecule  —  l*145^nj  Gj  +  Vnn  Vee 


(12) 


Eqs.  (11)  and  (12)  can  be  regarded  as  approximate  Kohn-Sham  analogues  of  the  exact 
Hartree-Fock  eqs.  (8)  and  (7).  Using  the  data  in  Tables  1  and  4,  eq.  (1 1)  reproduces  the 
respective  B3  total  energies  to  a  notable  degree  of  accuracy,  the  average  absolute  deviation  being 

only  0.20%. 


5.  Summary 

The  key  points  to  come  out  of  this  work  are  the  following: 

( 1 )  Eq.  (6)  predicts  Kohn-Sham  molecular  energies  to  within  an  average  0.70%  of  the  actual 

values. 

(2)  Eq.  (10),  with  appropriate  values  of  a,  predicts  Kohn-Sham  molecular  energies  to  within  an 
average  1.6%. 

(3)  Eqs.  (11)  and  (12)  are  approximate  B3  Kohn-Sham  analogues  of  the  exact  Hartree-Fock  eqs. 
(8)  and  (7).  Eq.  (11)  predicts  B3  Kohn-Sham  molecular  energies  to  within  an  average 

0.20%. 
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Figure  Captions 


Figure  1.  Relationship  between  Emoi/^T  n;  ej  and  percent  deviation  from  eq.  (6),  for 
HF/6-31+G**  results.  The  correlation  coefficient  is  0.9996. 

Figure  2.  Relationship  between  Emoi  Uj  ^nd  percent  deviation  from  eq.  (6),  for 
BLYP/6-31+G**  results.  The  correlation  coefficient  is  0.979. 

Figure  3.  Relationship  between  Emoi/^  nj  £j  and  percent  deviation  from  eq.  (6),  for 
B3P86/6-31 1++G(3df,3pd)  results.  The  correlation  coefficient  is  0.969. 


Table  1.  Energy  quantities  (in  hartrees)  computed  by  various  Hartree-Fock  and  Kohn-Sham  procedures. 
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Table  2.  Deviations  of  computed  results  from  eq.  (6):  %  difference  =  Emoj  -  —  (Vne  +  2V 
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Table  3.  Computed  values  of  mol 
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Table  4.  Computed  energy  quantities  (in  hartrees)  for  benzene,  aniline  and  nitrobenzene. 


Method 


Molecule 

Quantity 

HF/6-31+G** 

B3P86/6-31+G** 

c6h6 

Xni  €i 

-155.21902 

-138.03834 

Vne 

-945.88064 

-944.02137 

Vnn 

204.58471 

203.47523 

Emol 

-230.72183 

-233.02088 

^(Vne  +2Vnn) 

-230.01909 

-230.17325 

%  difference3 

-0.30% 

-1.22% 

Emol/21ni  ei 

1.486 

1.688 

c6h5nh2 

Xni  6i 

-190.97150 

-170.21654 

Vne 

-1211.67191 

-1209.36796 

V„n 

272.86409 

271.49366 

Emol 

-285.75786 

-288.52930 

|(V„,+2Vrai) 

-285.40446 

-285.59170 

%  difference3 

-0.12% 

-1.02% 

®mol/^ni 

1.496 

1.695 

c6h5no2 

Xni  €i 

-285.30524 

-256.28196 

Vne 

-1850.07489 

-1843.22742 

Vnn 

416.77305 

412.92240 

Emol 

-434.19709 

-437.91797 

|(V„e+2Vnn) 

-435.65520 

-436.02112 

%  difference3 

0.34% 

-0.43% 

^mol/Sni  ei 

1.522 

1.709 

aDefined  in  Table  2. 


Figure  Captions 


Figure  1.  Relationship  between  Emoj/^  n  j  Gj  and  percent  deviation  from  eq.  (6),  for 
HF/6-31+G**  results.  The  correlation  coefficient  is  0.9996. 

Figure  2.  Relationship  between  Emn|  /^,  n  j  Gj  and  percent  deviation  from  eq.  (6),  for 
BLYP/6-31+G**  results.  The  correlation  coefficient  is  0.979. 

Figure  3.  Relationship  between  Emol/ n,  Gj  and  percent  deviation  from  eq.  (6),  for 
B3P86/6-3 1 1++G(3df,3pd)  results.  The  correlation  coefficient  is  0.969. 
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